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ABSTRACT 

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted 
hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and 
Ha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the 
amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis 
of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for 
diffuse X-ray emission outside the optical disk, with a 3cr upper limit to the halo X-ray luminosity of 
4x 10 39 erg s _1 . An identical study of the less massive NGC 5170 also fails to detect any extraplanar X- 
ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical 
simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. 
For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities 
which are lower by a factor of ~ 2 compared to our earlier results reported by Toft et al. (2002). 
We find the new simulation predictions to be consistent with our observational constraints for both 
NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively 
star-forming spirals is in general probably associated with stellar activity in the disk. Observational 
results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints 
on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy 
formation and evolution. 

Subject headings: galaxies: formation — galaxies: halos — galaxies: individual (NGC 5170, NGC 5746) 
— galaxies: spiral — X-rays: galaxies 
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1. INTRODUCTION 

Estimates of the cosmic baryon fraction, defined as 
the ratio of baryonic to total mass in the Universe, 
can be combined with constraints on the integrated 
mass function of galaxies to infer that most baryons in 
the Un iverse are in a hot , diffuse form at the present 
epoch ([Balogh et al.l |2001[ ) . Cosmological simulations 
suggest that 30-40 per cent of all baryons reside in inter- 
galactic filaments of shock-heated gas with temperatures 
10 5 < T < 10 7 K, th e so-called warm /hot i ntergalactic 
medium, WHIM (e.g. ICen fc Ostrikei1ll999t iDave et all 
I200I . X -ray absorption studies have provided obser- 
vational evidence for such a componen t along the line of 
sight towards a number of quasars (e.g. Tripp et aL1l2000h 



iNicastro etai1l2005t [Savage et alj|20b5f )7 

While most of the WHIM baryons are predicted to re- 
side in structures of low overdensity, outside the dark 
mat ter halos of indiv idual galaxies and groups of galax- 
ies (jDave et al.ll200ll ). a potential repository for some of 
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the "hidden" baryons in the Universe could be extended 
halos of hot (~ 10 6 — 10 7 K) gas around individual g alax- 
ies, including spirals (e.g. lFukugita fe Pee bles 2006). For 
massive galaxies, cosmological simulations suggest that 
the total mass of these "external" galactic baryons is 
comparable to that of stars and col d gas in the galax- 
ies themselves (jSommer-Lar scn 2006) . In a cosmological 
context, observational support for such a scenario comes 
from the angular correlations between the galaxy distri- 
bution and the soft X-ray background, which indicate the 
presence of soft X-r ay emission f rom WHIM surrounding 
individual galaxies (Soltan 200(|). 

The presence of extended hot gaseous halos around op- 
tically bright elliptical galaxies is well establ ished from 
X-ray observations (e.g. lO'Sullivan et al.l l2001). The idea 
that disk galaxies could also be embedded in such halos 
is integral to man y semi-analytic a l mod e ls of disk galaxy 
forma tion (e.g. IWhite fc Reed 119781: IWhite fc Frenk 
19911 : ICole etaLl l2000t lHatton et al.l 120031 : iBower et al. 
2006 : ICroton et al.ll2006f K In these models, galaxy dark 



matter halos are assumed to grow as predicted by spher- 
ical infall models, with gas accreting continuously along 
with the dark matter. During infall into the dark mat- 
ter potential, the gas is heated, potentially to the halo 
virial temperature, subsequently cooling radiatively. If 
cooling is rapid, as is the case for characteristic gas 
temperatures T < 10 6 K and hence relatively shal- 
low potential wells, no accretion shock develops outside 
the evolving galactic disk. In these models, present- 
day spirals of total mass M < 10 11 M Q accrete gas 
which is predominantly in a cold, non-X-ray emitting 
phase at temperatures much lower than the virial tern- 
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perature of their halo (|Binnevl 
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2003 [lBinnevll2004tlKeres et all 
2006). This 'cold accretion' mode would be particu- 
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larly pronounced at high redshift, and would imply that 
most of the halo radiation is emitted as Lya emission 
close to the disk. Indeed, observational support for this 
scen ario is accumulating at both low and high redshift 



e.g., van der Hulst fe S ancisi 2004; Fraternali fc B innc 
2006: Nilsson et al. 2006: Sancisi et al. 2008: Smith et al 
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At higher galaxy masses and gas temperatures, the gas 
would emit at X-ray wavelengths but would cool less ef- 
ficiently, and is then assumed in the models to flow in 
more slowly than in the fast-cooling regime. A hot X- 
ray emitting halo can develop, in which the infalling gas 
is eventually deposited where dictated by its angular mo- 
mentum. As in the 'cold accretion' mode, the gas may 
ultimately condense to form stars in a rotationally sup- 
ported disk. In practice, thermal instabilities could act to 
locally reduce the cooling time of hot halo gas, so in ad- 
dition to a classical central cooling flow, cooling may also 
proceed via fragmentation of halo gas into smaller clouds 
throughout the halo. The net result may be the forma- 
tion of a multi-phase halo containing warm (T ~ 10 4 K), 
pressure-supported clouds embedded in a large-scale hot 
corona. These clouds may then rain down onto the 
disk, possibly akin to some of th e high- velocity Hi clouds 
seen around the Galaxy (e.g.. [Mailer fc Bullock! [20041: 
Kaufmann et al.ll200 6: Sommcr-Larscn 2006; Peck et all 
20081 : but see also iFraternali fc Binneyl l2008f) . In the 
scenario proposed by iMaller fc Bullock! (|2004f ). conden- 
sation and infall of such clouds would be very impor- 
tant for fueling MW-sized galaxy disks, although at least 
some numerical work suggests th at this mode of accre- 
tion is subdominant at z = (ISommer- Larsenl |2006; 
iKaufmann et alJl200d iPeek et al.ll2008l ). 

The general picture of massive disk galaxies embed- 
ded in hot gaseous halos has been backed by cosmologi- 
cal hydrodynam i c simulations of ga laxy formation (e.g., 
iToft et alJ [20021 : iKeres et al.ll2005h . thus lending sup- 
port to the general validity of the simplifying assump- 
tions underlying semi-analytic al models. For example, 
the simulations presented by IToft et all (|2002t ) showed 
the X-ray luminosity Lx of hot halo gas to scale strongly 
with disk circular velocity v c (defined here as the ro- 
tation velocity at 2.2 times th e disk scale-length; see 
ISommer-Larsen fc Dolgovll2001l for details), with Lx oc 
v^~ 7 , as also expected from semi-analytical models of 
galaxy formation. It is rele vant to note, howe ver, that all 
the massive galaxies in the lToft et all (|2002f ) study were 
formed in the old "standard" (O ro = 1, J7a = 0) cold dark 
matter cosmology, with a fixed primordial metal abun- 
dance and an assumed cosmic baryon fraction well be- 
low the currently accepted value. Thus, it seems timely 
to explore the hot halo predictions of more sophisticated 
higher-resolution simulations assuming the currently pre- 
vailing ACDM cosmology. This is one of the goals of this 
paper. 

The existing theoretical results clearly suggest that X- 
ray observations could provide a useful means of mapping 
any material accreting onto massive spirals at the present 
epoch. Direct observational verification of the above sce- 
nario has remained elusive, however. In a first study 
aimed at detecting cosmologically accreted gas halos 



around spirals, iBenson et al.1 (|2000[ ) compared ROSAT 
X-ray observations of three massive, highly inclined spi- 
rals to predictions of simple cooling flow models. The 
failure of these authors to detect any hot gas around 
the galaxies translated into upper limits on the halo X- 
ray luminosity which, for their best-constrained case of 
NGC2841, was found to be a factor of 30 below their 
model prediction. While their models were deliberately 
simplified versions of general semi-analytical ones, the re- 
sult persists that some hydrodynamical simulations and 
simple a nalytical models overpredict actual halo lu mi- 
nosities (jGovernato et al.l l2004t IBenson et aI1l2000l and 
references therein), though the exact magnitude of this 
discrepancy has yet to be established. 

Extended X-ray halos have been detected around 
numerous lat e -type star-forming s p irals (see 
IStrickland et alJ l2004at ITullmann etai] l2006al and 
references therein). However, these galaxies typically 
have at least moderately high star formation rates 
(> 1 M yr" 1 ), with the detected halo X -ray emission 
coinciding with extraplanar Ha and radio continuum 
emission. The standard interpretation is that these 
multi-phase halos have been generated by star formation 
activity in the disk, with supernovae expelling hot, X-ray 
emitting gas, warm and/or photoionized Ha-emitting 
material, and radio syn chrotron emitting cosmic ray 
electr ons from the disk (IWang et al.l 120011. [2001 lEhlei 
12004 IStrickland et all l2004aUbl : ITullmann et al.ll2006aft . 
Unlike the expectation for gravitationally heated gas, 
the X-ray luminosity of the gas in the halos of these 
galaxies is found to correlate with various measures of 
the disk star formation rate and with the energy input 
rate by supernovae (SN), but — even when including 
the diffuse X-ray luminosity of the disk — not with 
baryoni c or Hi mass in the disk (jTiillmann et al.ll2006bl: 
see also lSun et al.ll2007l ). 

Diffuse X-ray emission has also been detected around 
massive early-type (Sa-S b) spirals such as NGC 4594 
(the "Sombrero" galaxy) ([Li et al.ll200l I2007D . In these 
cases, the emission is not obviously associated with star 
formation in the disk, but the total diffuse X-ray lumi- 
nosity from the disk and halo still falls well below expec- 
tations from cosm ological accretion on the basis of the 
IToft et all (|2002ft results. It has been speculated that 
much of the detected extraplanar emission in these galax- 
ies is heated by S N la activity in the disk rather than by 
gravity (see, e.g.,[Wang 2003). 

Summarizing, the hot halos observed so far around 
spirals seem to result mainly as a consequence of stel- 
lar feedback in the disk rather than infalling intergalac- 
tic material. A possible exception is the Milky Way it- 
self, a moderately star-forming galaxy for which indi- 
rect arguments can be used to ascertain the presence of 
an extended rese rvoir of hot gas possibly related to a 
gaseo u s halo (e.g. iMoore fc Davisl[l9 94: Quili s fc Moore 
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Scmbac h et al.l l2003t iBregman fc Llovd-Davies 
IWangi l2007f T Despite this, it remains unclear 
whether any of this putative halo gas is actually cool- 
ing out to accrete onto the Milky Way disk, and if so, 
to what extent it can replenish the gas currently being 
consumed in the disk by star formation. 

Sensitive X-ray observations of carefully selected sys- 
tems are required to settle the question of whether sub- 
stantial hot gas halos associated with infalling material 
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exist around quiescent spirals at all. This would pro- 
vide an important consistency check on the assumptions 
underlying existing disk galaxy formation models. Mo- 
tivated by this, we initiated a Chandra search for hot 
X-ray gas surrounding massive, quiescent disk galax- 
ies. As part of this program, the discovery of hot halo 
gas around NGC 5746, with a total X -ray luminosity 
of ~ 4 x 10 39 erg s _1 , was reported by iPedersen et all 
(2006). In the present paper, we re- visit this issue in 
more detail using updated Chandra calibration data, de- 
scribe our corresponding Chandra results for the less 
massive galaxy NGC 5170, and present new results from 
cosmological simulations of galaxy formation and evolu- 
tion, for direct comparison to the obse rvations. Contrary 
to the claim in IPedersen et alj (2006) , significant extra- 
planar emission surrounding NGC 5746 is not detected 
in our updated analysis. 

In § [2] we outline our target selection and the analy- 
sis of X-ray and optical data. Results for both galaxies 
are presented in § [3] and compared to predictions of cos- 
mological simulations and to results for generally more 
actively star-forming spirals in § [4] Some implications 
of our observational and theoretical results are discussed 
in § [SJ and our findings and conclusions are summarized 
in § H A Hubble constant of H = 73 km s" 1 Mpc -1 
is assumed throughout. The distances to NGC 5746 and 
NGC 5170 as listed in NED are then 26.8 and 24.8 Mpc, 
respectively, with 1 arcmin corresponding to ~ 7.7 and 
~ 7.1 kpc. Unless otherwise stated, uncertainties are 
reported at the 68 per cent confidence level. 

2. GALAXY SAMPLE AND DATA ANALYSIS 
2.1. Target Selection 

The primary objective of this study is to map and char- 
acterize any diffuse X-ray gas surrounding quiescent disk 
galaxies. Our starting point w as therefore the simula- 
tion results of iToft et al.1 ([2002), in particular the result 
that halo X-ray luminosity Lx should scale strongly with 
disk circular velocity v c . Hence, our target sample was 
drawn by searching the NED and HyperLeda databases 
for nearby spiral galaxies with circula r velocity y c > 28 
km s -1 . According to the results of IToft et all (|2002ft . 
this should ensure that the target galaxies are massive 
enough to generate and retain a large-scale X-ray halo 
detectable within reasonable Chandra exposure times. A 
further advantage of large v c lies in the fact that the 
efficiency with which SN outflows can escape the disk 
and complicate the interpretation of our results may be 
expected to decrease with increasing disk mass. 

In addition, we required our targets to meet the fol- 
lowing criteria: i) inclination i > 80°, to enable a clean 
separation between disk and halo gas; ii) distance D < 
40 Mpc, for similar reasons, and to obtain a sufficiently 
strong X-ray signal; iii) Galactic latitude |6| > 30°, 
to have a low column density of absorbing foreground 
Galactic hydrogen, given that any halo emission is ex- 
pected to be fairly soft and thus subject to significant ab- 
sorption by intervening material (thi s criterion excl uded 
a galaxy like NGC 2613, studied bv lLi et al.|[2006l ); iv) 
the galaxies should not be a member of a known X-ray 
bright group or cluster, nor show any significant signs 
of interaction, so as to include only undisturbed galaxies 
and avoid significant contamination from hot gas residing 



TABLE 1 

Salient parameters of the observed disk galaxies 



Parameter 


NGC 5746 


NGC 5170 


D (Mpc) 


26.8 


24.8 


RA (J2000) 


14' l 44 m 56?00 


13 h 29 m 48!83 


Dec (J2000) 


+01°57'17'.'l 


-17°57'59'.'4 


Hubble type 


SBb 


Sc 


v c (km s ) 


318 ± 10 


247 ± 3 


D25 (arcmin) 


(6.92, 1.20) 


(8.32, 1.20) 


M B 


-21.79 


-21.18 


M K 


-25.26 


-24.35 


£tr (L©) 


6.2 x 10 9 


3.7 x 10 9 


SFR (M s yr" 1 ) 


0.9 ±0.2 


0.5 ±0.1 


Inclination 


84° 


90° 


N H (cm- 2 ) 


3.3 x 10 20 


6.9 x 10 20 


Chandra obs. ID 


3929 


3928 


Chandra obs. date 


2003-04-11 


2003-05-18 


Chandra exp. (ks) 


35.7 


31.9 



Note. — Distances, positions, and magnitudes taken from the 
NASA/IPAC Extragalactic Database (NED). Disk circular velocity 
"ci £*25 (major and minor diameter of the ellipse outlining a B— 
band isophotal level of 25 mag arcsec - 2 ), and inclination from the 
HyperLeda database. Ljr is the total (8-1000 )im) infrared lumi- 
nosity , based on IR A S fluxes and the relation of [Sanders & Mirabell 
(1996), and SFR the estimated star formation rates (see text for 
details). Absorbing colum n density JVh is the Galactic value from 
IDickev fc Lockmanl Jl990T) . 



in an ambient intragroup/-cluster medium; v) in order 
to also minimize contamination from outflows of hot gas, 
the galaxies should be relatively quiescent, i.e. not show 
evidence for significant activity related to a starburst or 
an active galactic nucleus (AGN). This excluded galax- 
ies such as NGC 891, which displays extraplanar X-ray, 
Ha, and rad io emission associat ed with stellar feedback 
activity (see lTemple et aLll2005l and references therein), 
and the strongly bulge-dominated Seyfert NGC 4594 (the 
' Somb rero'), which hosts a ~ 10 9 M w low-lum inosity 
AGN dFabbiano fc Juddfl997l iBendo et al.ll2006h . 

Imposing these strict criteria, only one galaxy re- 
mained, NGC 5746. We therefore relaxed the require- 
ment on v c to include a second target, NGC 5170, which, 
although still massive, could act as a consistency check 
on our analysis procedure. NGC 5746 is an SBb spiral 
with a circular velocity v c — 318 ± 10 km s _1 (Hyper- 
Leda database), and NGC 5170 an S c spiral with a signif - 
icantly lower v c of 247 ± 3 km s _1 ([Kregel et al.ll2004h . 
We note that our selection was not restricted to par- 
ticular Hubble types, but given the v c > 280 km s _1 
criterion, the selection method is likely to bias against 
late-type spirals. Since such galaxies are likely to display 
higher star formation activity at present than early-type 
ones, this criterion is in turn likely to favor relatively 
quiescent spirals, as desired. Salient parameters of the 
observed galaxies are listed in Table [TJ 

As the goal of this study is to constrain the amount 
of infalling hot halo gas, the requirement of low lev- 
els of disk activity is particularly important for ensur- 
ing minimal contamination of the extraplanar X-ray flux 
from outflowing hot gas. In terms of their star forma- 
tion rate (SFR), both galaxies are indeed fairly quies- 
cent, as already evidenced by the fact that none of them 
is in cluded in the revise d catalog of bright IRAS galax- 
ies ([Sanders et al.1 12003T ), despite their close proximity 
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and large mass. The ir IRAS 60 and 100 fim fluxes S^o 
and ^nn, taken from [Moslrkelal] (fl990l) for NGC 5746 
and lRice et ail (| 19881) for NGC 5170, translate into global 
SFRs of 0.9 ± 0.2 (NGC 5746) and 0.5 ± 0.1 yr" 1 
(NGC 5170 ], using the rela tions of ISanders fe Mirabel 
(1996) and iKennicuttl ((1998). This is equivalent to spe- 
cific (area-normalized) SFRs of ~ 4 x 10~ 4 (NGC 5746) 
and - 2 x 10~ 4 M yr" 1 kpc~ 2 (NGC 5170), assum- 
ing a circular disk of radius equal to the semi-major 
axis of D25- For comparison, the moderately quies- 
cent Mil ky- Way disk displays an SFR of ~ 2 M Q yr _1 
(see e.g. iCasuso fc Beckmanl 120041 ). corresponding to ~ 
3 x 10 -3 Mq yr _1 kpc~ 2 . Both galaxies furthermore dis- 
play low far-infrared 'temperatures' of Seo/Swo ~ 0.15 
(NGC 5746) and 0.22 (NGC 5170), and their far-infrared 
luminosities £fir imply " mass- normalized" star forma- 
tion rates Lfir/^b « 0.06, further indicative of lo w star 
formation activity (see e.g. lRead fc Ponm an 2001 for de- 
tails) . 

From a Chandra study of extraplanar X-ray emis- 
sion around nearby disk galaxies cov ering a range of 
morpho logies and star formation rates, IStrickland et al.l 
(2004b]) derived a criterion for gas blow-out by super- 
novae based on the observed properties of the inter- 
stellar medium in typical spirals. They find on the- 
oretical grounds that in order to have a collection of 
near-simultaneous SN expel gas from the disk, the re- 
quirement -Fsn,fir,d 25 > 25 SN Myr -1 kpc~ 2 must 
be satisfied. In their notation, -Fsn,fir,.D25 is the rate 
T^-SN = 0.2-Lir/IO 11 L q of core-collapse SN per disk 
"area" -D| 5 . Indeed, for the eight out of ten spirals in 
their sample with detected extraplanar X-ray emission, 
they find F S n,fir,d 25 > 40 SN Myr" 1 kpc~ 2 . Figure CD 
co mpares the core-c o llapse SN rates of the galaxies in 
the IStrickland et alj (|2004al lH) sample to those inferred 
for our target galaxie s using the same meth ods. The 
additional galaxies of Tiillman n et al.l (|2006al ) from a 
comparable XMM-Newton study have also been plotted. 
Both our targets clearly lie in a regime in which ejection 
of detectable hot gas by SN would not be anticipated. 
While these results neglect the contribution from type la 
SNe, low specific SN la rates of w 4.7 (NGC 5746) and 
3.3 SN Myr" 1 kpc~ 2 (NGC 5170 ) can be estimated, us - 
ing the prescription adopted bv iMannucci et all (|2005l ) 
and the K- and £>-band magnitudes listed in Table [TJ 
Adding these to the results in Fig. Q] clearly still places 
bo th galaxies well b elow t he blowout criterion derived 
bv IStrickland et all (|2004bl ). These results suggest that 
any extraplanar emission surrounding our targets will 
be largely uncontaminated by emission from supernova- 
driven outflows of hot gas, enabling reliable constraints 
on the physical properties of any accreting material. 

2.2. Observations and Analysis 

Both disk galaxies were observed by Chandra with the 
ACIS-I array as aimpoint and with the CCD's at a tem- 
perature of —120° C. Data were telemetered in Very 
Faint mode which allows for superior background sup- 
pression relative to standard Faint mode. To exploit 
this, the data were reprocessed and background screened 
using CIAO v3.4 with CALDB 3.4.1. Bad pixels were 
screened out using the bad pixel map provided by the 
reduction pipeline, and remaining events were grade fil- 
tered, excluding ASCA grades 1, 5, and 7. Periods of 
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Fig. 1. — Comparison of the specific core-collapse SN 
rates estimated for o ur target galaxie s (stars ) to those in the 
disk galaxy sample of Strickland ct al. ( 2004a. b) (triangles) and 
ITiillmann et all d2006al ) (squares) . Galaxies with detected extra- 
planar X-ray emission (in all cases attributed to stellar feedback) 
are represented by filled symbols. The dashed line marks the ex- 
pected critical SN rate, ~ 25 SN Myr -1 kpc -2 , above which SN 
are expectedly capable of expelling gas from the disk. 



high background on the ACIS-I chips were filtered us- 
ing 3cr clipping of the 0.3-12 keV lightcurves extracted 
in off-source regions in 259-s bins. Resulting lightcurves 
showed no strong flaring periods, leaving a total of 35.7 ks 
(NGC 5746) and 31.9 ks (NGC 5170) of cleaned exposure 
time. 

Point source searches were carried out with the CIAO 
task 'wavdetect' using a range of scales and detection 
thresholds, and results were combined. A total of 126 
(NGC 5746) and 111 (NGC 5170) point sources were de- 
tected in the two observations, of which 20 and 17 
sources, respectively, are located inside the D25 ellipses. 
Point source extents were quantified using the 4cr de- 
tection ellipses from 'wavdetect', and these regions were 
masked out in the analysis of diffuse emission. Further 
details of the X-ray analysis are described below. 

In addition to the X-ray data, we also obtained Ha im- 
ages of both galaxies with the Danish 1.54-m telescope 
at La Silla, Chile, for total exposures of 120 (NGC 5746) 
and 160 min (NGC 5170). In both cases, 20 min R- 
band images were also taken for subsequent continuum- 
subtraction. The frames were bias-subtracted, flat- 
fielded, and median-combined using standard iraf pro- 
cedures. Using the i?-band image of each galaxy, we 
estimated and subtracted the continuum emission con- 
tribution to the Ha image as follows. The two combined 
images were normalized to the same exposure time, and 
the i?-band image was smoothed with a Gaussian ker- 
nel to match the broader point spread function of the 
Ha image. The smoothed i?-band image was then scaled 
to match the narrower width of the Ha filter and sub- 
tracted from the Ha image. For NGC 5746, the back- 
ground in the resulting continuum-subtracted Ha image 
had a large-scale gradient caused by stray light from a 
star just outside the field of vie w. We modeled this us- 
ing the SExtractor software ([Bertin fc Arnoutslll996f) 
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Fig. 2. — 0.3-2 keV exposure-corrected X-ray contours overlayed on our Ha images of both galaxies. Contours have been smoothed with 
a Gaussian of width cr = 20", and are logarithmically spaced over a decade beginning at 1.5 times the background level. The D25 ellipse 
is outlined by a dashed contour in each case. 



with options set to save a full resolution interpolated 
background map. Parameters were optimized to gener- 
ate a map sufficiently fine to represent large scale vari- 
ations, and sufficiently coarse so as not to include local 
non-background structures (such as an Ha halo around 
the galaxy). The resulting background map was then 
subtracted from the Ha image. 

3. RESULTS 

3.f. X-ray and Ha Emission 

In order to aid the search for diffuse X-ray emission 
around either galaxy, smoothed, exposure-corrected 0.3- 
2 keV images were created for both using two separate 
methods. One involved smoothing the images with a 
simple Gaussian kernel, while in the other we generated 
backgroun d-subtracted imag e s foll owing the procedure 
outlined in iRasmussen et all (|2006f h In either case, the 
resulting images were exposure-corrected using similarly 
smoothed exposure maps. The overall impression con- 
veyed by these images is qualitatively similar, suggesting 
very little diffuse emission beyond the optical disk of ei- 
ther galaxy. Given the expected weakness of any halo sig- 
nal, we here only discuss the simple Gaussian smoothed 
images, so as to accommodate any concerns that adap- 
tive smoothing may introduce spurious features for these 
low-S/N data. Fig. [2] shows our Ha images of the cen- 
tral 12 x 12 arcmin region around each galaxy, with the 
Gaussian smoothed X-ray contours overlayed, and with 
the optical extent of the galaxies outlined by their re- 
spective D25 ellipses. Clearly, both diffuse and point-like 
X-ray emission is detected within the optical disks of 
both galaxies, but there is no immediate indication of 
any diffuse X-ray or Ha emission extending well beyond 
the disks. 

In order to perform a more quantitative search for dif- 
fuse X-ray emission around the galaxies, and confirm 
the overall impression conveyed by Fig. [21 we generated 
0.3-2 keV surface brightness profiles of the unsmoothed, 
exposure-corrected emission. Point sources, chip gaps, 



and emission inside D25 were masked out. Radial profiles 
extracted within r = 40 kpc (r rj 5 arcmin for NGC 5746) 
are shown in Figure [3l binned into ten equal-sized radial 
bins for improved S/N, and with the background level in 
the data evaluated from a surrounding (r = 5.5-7.5 ar- 
cmin) point source-excised annulus. 

These profiles support the conclusion that no signifi- 
cant detectable X-ray emission surrounds either galaxy 
on large scales. We confirmed that this result does not 
depend on the specific choice of radial binning scheme. 
While the NGC 5170 profile shows emission outside D25 
which is everywhere consistent with a constant level 
equal to the estimated mean background level, the profile 
for NGC 5746 does show some deviations from a uniform 
level. Excess emission is detected in the radial bin im- 
mediately outside the optical disk at 2.2cr significance, 
suggesting the presence of faint extraplanar emission at 
~ 5-10 kpc from the galactic center. We investigate this 
in more detail in the inset in Fig. [3^i, which displays the 
NGC 5746 profile extracted in bins perpendicular to the 
disk. This plot confirms that the excess is not an artifact 
of the adopted binning scheme, while also indicating that 
the excess may be largely confined to the eastern side of 
the disk. However, since the signal is only nominally sig- 
nificant at the ~ 2cr level and is comprised of only 15 net 
photons, we do not consider it robust evidence of diffuse 
halo emission associated with NGC 5746. It is also worth 
remarking that a comparably significant deficit of emis- 
sion is seen at larger radii (r s=s 3.5-4 arcmin), suggest- 
ing that the excess may be a spurious feature introduced 
by local fluctuations in the background level, perhaps 
caused by the presence of point sources just below our 
detection limit. 

Inside the circular aperture underlying Fig. [3) i.e. with 
the D25 ellipse excluded, we detect a total of 38 net pho- 
tons in 0.3-2 keV for NGC 5746, corresponding to emis- 
sion at just 0.8cr above the background. The correspond- 
ing result for NGC 5170 shows a deficit of 101 photons 
(significant at 2.1a) relative to the surrounding back- 
ground. Qualitatively similar results are obtained in 
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Fig. 3. — 0.3-2 keV exposure-corrected radial surface brightness 
profiles of (a) NGC 5746 and (b) NGC 5170. Dashed vertical lines 
outline the maximum extent of the D25 ellipse minor axis, and 
shaded regions mark the 90% confidence uncertainties on the local 
background level. Point sources, chip gaps, and emission within 
the D25 ellipse itself have been masked out. Inset in (a) shows 
the NGC 5746 profile (same units, with east to the left) extracted 
perpendicular to the disk, with D25 again masked out. 



narrower energy bands such as 0.5-1.5 keV, in which 
the number of source photons is obviously lower but the 
ACIS-I S/N from any absorbed T ~ 0.2-0.3 keV plasma 
should be higher. Hence, the immediate conclusion is 
that no statistically significant diffuse X-ray emission is 
detected outside the optical disk of either galaxy on these 
scales. 

In other cases where diffuse X-ray emission has been 
reported around spirals, this emission is typically ac- 
companied by extended Ha an d radio halos, with sim- 
ilar overall morpholog ies (e.g., IStrickland et al.1 [2004a; 
Tiillmann et al. 2006a). Although there is no indication 
in Fig. [2] of significant extraplanar Ha emission around 
NGC 5746, we investigate this possibility in more detail 
in Fig. [31 where we show the surface brightness profile 
of the Ha emission perpendicular to the disk of this 
galaxy. As can be seen, there is no evidence for signifi- 
cant amounts of Ha gas outside the optical extent of the 
disk. Note that the dip in the Ha profile 10 arcsec east of 
the D25 center is due to a dust lane also visible in broad- 
band optical imag es. For NGC 5 170, a simi lar study has 
been performed bv lRossa &: Dettmarl (|2000T ). who find no 
evidence for extraplanar diffuse Ha emission surrounding 



this galaxy in an Ha exposure of similar depth to ours. 
There is also no indication of significant extraplanar ra- 
dio emission around either galaxy in the public NRAO 
VLA Sky Survey (NVSS) 1.4 GHz data of iCondon etafl 
(1998). For NGC 5746, this conclusion is further sup- 
ported by the NVSS flux density profile extracted from 
these data, as shown in Fig. [3b. These results provide 
further testimony to the low star formation activity of 
both our galaxies, and support the assumption based on 
Fig. [T] that our constraints on X-ray emission associated 
with any gas accreting onto these galaxies are not subject 
to serious contamination by emission from disk outflows. 

3.2. Differences with Respect to \Pedersen et al\ \2006i) 

The absence of significant detectable diffuse emission 
around NGC 5 7 46 co ntrasts with the earlier result of 
iPedersen et al.l (|2006l ) in which a 4tr detection of halo 
emission was claimed, with a 0.3-2 keV luminosity Lx = 
10 39 erg s _1 and a plasma temperature of 
a20 keV. Although the modest significance 



T 



0.56] 



of the detection renders this result tentative rather than 
conclusive, the difference with respect to our present con- 
clusion demands an attempt at an explanation. We spec- 
ulate that this discrepancy arises due to a change in the 
calibration of the time-dependent degradation of the low- 
energy ACIS CCD quan tum efficiency, a possibility also 
advocated by I Wang! (|2007h . This degradation is believed 
to be caused by a build-up of molecular ice on the Chan- 
dra/ ACIS optical blocking filters or CCD's 1 , an effect 
which was known but still no t well understo o d at the 
time of our original analysis in IPedersen et a l. (2006). 

The properties of the contaminant has remained the 
focus of continued measurements and is considerably bet- 
ter understood at present. In particular, a model of the 
spatial dependence of the contaminant was introduced 
with the calib ration data in CALD B v. 3.0. However, 
the analysis of IPedersen et a l. (2006) was performed us- 
ing CALDB v. 2.26, with which the contamination could 
only be appropriately corrected for in the case of an on- 
axis source. As contamination is observed to increase 
towards the edges of the ACIS array, not correctly ac- 
counting for this spatial variation will introduce a spuri- 
ous signal in exposure-corrected low-energy ACIS-I im- 
ages, appearing as faint diffuse emission whose surface 
brightness decreases with distance from the CCD aim- 
point. It will also lead one to underestimate the soft 
X-ray background at large off-axis angles. 

We conducted a few tests to explore whether this 
effect could account for the halo signal detected by 
IPedersen et al.l (|2006h . At a photon energy E — 0.7 keV, 
close to the flux peak energy of a T w .6 keV thermal 
plasm a (the halo temperature inferred bv IPedersen et all 
l2006h . the additional optical depth towards the edges of 
the CCD arrays can be 0.3-0.5 for data taken around 
the time of our observations (see previous footnote), im- 
plying an apparent reduction in surface brightness at 
large off-axis angles to 60-75% of the on-axis value. The 
best-fit (halo plus background) surf ace brightness model 
reported by IPedersen et al.l (|2006h in fact declines at 
large distances from the NGC 5746 disk to about 75% 
of the value seen immediately outside the optical disk 

1 http:/ /cxc. harvard.edu/cal/Acis/Cal_prods/qcDeg/ 
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Fig. 4. — (a) Continuum-subtracted Ha surface brightness profile of NGC 5746, extracted perpendicular to the disk within the same 
aperture as that underlying the inset in Fig. [3^. (again with east to the left). Bright stars have been masked out. (b) Corresponding mean 
flux density of 1.4 GHz continuum emission within the same aperture, in bins corresponding to the VLA beam size of 45 arcsec (FWHM). 
Bright point sources have been masked out. Uncertainties represent the standard deviation within each bin, and dotted lines outline the 
la errors on the noise level. Dashed vertical lines are as in Fig. \3\ 



(at \z\ ~ 1.5'), consistent with the expectation for the 
spatial variation of the contaminant. 

As a further test of this hyp othesis, we re-consider the 
spectral analysis described by iPedersen et al.1 (|2006D , in 
which a double-subtraction approach was employed for 
the extraction of a background spectrum. This involves 
extracting a background spectrum for the same detector 
region as the source region, using blank-sky background 
data prepared in a similar manner to the source data. 
This procedure is then repeated for a large-radius an- 
nulus, in order to obtain a residual background spec- 
trum that helps to account for any difference in soft 
X-ray background between source- and blank-sky data. 
This residual spectrum is then added to the source-region 
blank-sky spectrum to produce a composite background 
sp ectrum. We her e repea ted the exact same procedure as 
in IPedersen et al.l (|2006f ) , while employing the more re- 
cent calibration files included with CALDB 3.4.1. In or- 
der to mimic the effect of underestimating the soft X-ray 
background at large radii, as implied by our hypothesis, 
the normalization of the soft residual spectrum was artifi- 
cially scaled down so as to ac hieve ~ 200 net c ounts in the 
source region as reported by IPedersen et ail ((2006) . The 
resulting background-subtracted source spectrum was 
then fitted precisely as in our original analysis, yielding a 
best-fit temperature and unabsorbed 0.3-2 keV luminos- 
ity of T = 0.561^5 keV and L x = 4.0±J;| x 10 39 erg s -1 , 
respectively, with y 2 — 7.1 for 9 degrees of freedom. 
These values are str ikingly similar to those obtained by 
IPedersen etHI (|2006l) . 

In summary, both the spatial and spectral properties 
of the diffuse X-ray h alo around NGC 5746 reported by 
IPedersen et al.l (|2006f ) are consistent with the hypothesis 
that the apparent halo signal was an instrumental arti- 
fact caused by the spatial variation in the degradation of 
the low-energy ACIS quantum efficiency, an effect which 
our original analysis could not take into account. 

The question then arises as to why this calibration 
problem did not lead to a similar spurious detection of a 
hot halo around NGC 5170, given that the two galaxies 
were observed with ACIS-I within one month of each 



other and with comparable exposure times. The low 
signal-to-noise ratio outside the disk in both data sets, 
including that of the original halo "detection" around 
NGC 5746, precludes robust conclusions in this regard 
however, so we can only offer some general considera- 
tions. We note that we find no statistically significant 
differences in the background level between the two data 
sets in any of a dozen different energy bands or large- 
radius background regions considered, so global back- 
ground differences are unlikely to play a major role. 
The slightly different data preparation criteria applied 
in our updated analysis could have some impact though, 
changes in the calibration data aside. For example, our 
lightcurve screening employs somewhat larger time bins, 
and has removed a previously undetected minor flare 
at the end of the NGC 5170 observation. Our current 
point source search is also slightly more sophisticated, 
and has resulted in the identification of ten additional 
point-like sources in the NGC 5170 data (but none such 
for NGC 5746), of which seven are outside the disk. In 
combination with the slightly shorter exposure time ob- 
tained for NGC 5170, these issues may have played a role 
in suppressing any faint spurious signal that could mimic 
diffuse X-ray emission around this galaxy. In addition, 
we cannot exclude the possibility that the spurious halo 
signal around NGC 5746 had an additional contribution 
not present in the NGC 5170 data, e.g. from faint un- 
resolved sources, such as the one hinted at immediately 
cast of the NGC 5746 disk in Figs. 1 and Ek 

3.3. Constraints on Halo Gas Properties 

Based on the absence of detectable extraplanar emis- 
sion around either of our targets, we can constrain some 
basic properties of any hot halo gas surrounding these 
galaxies. In this context, it is important not to confuse 
the signature of cooling halo gas with X-ray emission 
related to processes in the disk, such as the integrated 
emission from unresolved discrete sources and from hot 
gas heated by supernovae and stellar winds. The opti- 
cal disk of either galaxy is thus excluded in the following 
analysis. On the other hand, the apparent absence of 
significant X-ray, Ha, or radio emission outside the opti- 
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TABLE 2 

Derived constraints on hot halo gas 



Galaxy 


T 


Z 


L X 


(n e ) 


A^hot (*cool) 


(Jli) 




(keV) 


(Z ) 


(erg s _1 ) 


(cm" 3 ) 


(M ) (Gyr) 


Mq yr" 1 


NGC 5746 


0.3 (fixed) 


0.2 (fixed) 


< 4.0 X 10 3a 


< 3.5 x 10" 4 


< 2.3 x 10 9 > 7.9 


< 0.3 


NGC5170 


0.2 (fixed) 


0.2 (fixed) 


< 7.5 x 10 39 


< 5.3 x 10" 4 


< 3.4 x 10 9 > 2.3 


< 1.5 



Note. — All constraints are given at 3<r significance and are derived within r = 40 kpc, with the D25 ellipses excluded. Lx is given in 
the 0.3—2 keV band, corrected for Galactic absorption, and (M) M^ot / (t cool } is t ne mean cool-out rate of hot gas within the assumed 
aperture. 



cal disks, combined with the low specific star formation 
and supernova rates of the galaxies as discussed in § [3J 
suggest minimal contamination from outflows of hot gas 
driven beyond the optical disk of either galaxy by stellar 
feedback. Hence, we base our constraints on the amount 
of accreted hot halo gas on the net count rate limits eval- 
uated inside the circular apertures underlying Fig. [3] i.e. 
a region extending to r = 40 kpc from immediately out- 
side the D25 ellipse of each galaxy. 

In order to convert observed count rate limits into con- 
straints on the halo X-ray flux, we assume that the spec- 
trum of any hot halo gas can be described by a ther- 
mal plasma with a given temperature, metal abundance, 
and volume filling factor (adopting the APEC model in 
xspec vll.3). In the absence of direct observational in- 
formation on these quantities, we resort to more indirect 
observational constraints and to simulation results. For 
example, cosmological simulations predict a remarkably 
tight correlation between the temperature of i nfalling hot 
gas a nd v c of the disk (§ U below; see also iToft et all 
2002) , with T being in good agreement with a simple ex- 
pectation for the virial temperature T v - lr rj (l/2)/j,m p v^ 
of the underlying dark matter halo. We therefore nom- 
inally assume T = T v ; r w 0.3 and « 0.2 keV for any 
halo gas around NGC 5746 and NGC 5170, respectively, 
based on the disk circular velocities listed in Table Q] As 
further justification of this choice, we note that gas much 
colder than T V i r would mainly emit at UV wavelengths 
and likely be subject to rapid cool-out to T < 10 4 K, 
while gas much hotter would be subject to evaporation 
from the galactic gravitational potential. These situa- 
tions are manifestly distinct from the case of a gravita- 
tionally confined hot halo composed of infalling shock- 
heated material as seen in our simulations, but it is the 
amount and properties of any such material that is the 
focus of our attention. 

We further note that outside the galactic disks them- 
selves, the emission-weighted metal abundance of the X- 
ray gas in these simulations displays a mean and stan- 
dard deviation of Z — 0.2 ± 0.1 Z Q , which is consis- 
tent with observatio nal results for the lo w-redshift in- 
tergalactic medi um (iDanforth et al.ll2006l), fo r the out- 
skirts of groups asmussen fc Po nman 2007]) and clus- 
ters of galaxies fe.g.. iTamura et al.ll2004f) . and for the 
intergalactic fil aments from which the se structures ac- 
crete material (jFinoguenov et aL1l2003l ). Hence, for the 
purpose of our spectral analysis, we here consider a fidu- 
cial value of Z = 0.2 Z for the halo metal abundance 
(while noting that the mass-weighted abundance in the 
simulations is substantially lower, as discussed in § I4.3j) . 
Finally, for the volume filling factor 77 of hot g value 
of unity is assumed for simplicity, since our simulations 



(§S} indicate high values of r\ w 0.8-1. 

Using these assumptions and the appropriate ACIS-I 
response files, observed count rate limits were translated 
into corresponding constraints on the APEC plasma 
model normalization A, 

10~ 14 f 

A = 4^ J neHH dV > (1) 

where D is the assumed distance (Table [lj, and n e and 
71h are the number densities of electrons and hydrogen 
atoms, respectively. The integral represents the emission 
integral I c = J n e riH dV « Tfrv^V of the hot halo gas in- 
side the covered volume V . From the resulting constraint 
on I Cl we obtained upper limits to the mean hot gas den- 
sity (n e ) ~ (Ic/Vrj) 1 /' 2 w (I c /V) 1/2 and total hot gas 
mass Mhot ~ rn p {n e )V inside V . For V itself, a spherical 
volume of r — 40 kpc was assumed, with the D25 "cylin- 
der" excluded alo ng the line of sight. U s ing th e cooling 
curves A(T,Z) of lSutherland fc Dopital (pM) . we also 
derived constraints on the mean cooling time (t coo \) « 
3fcT/(A(n e }) and cool-out rate (M) w Mh t/(*cool) of 
hot gas inside this volume. 

Table [5] summarizes the derived constraints on halo 
properties for both galaxies, where all values have been 
obtained from the 3tr upper limits on the observed halo 
count rates using the nominal values of T, Z , and 77 dis- 
cussed above. Note that although the two galaxies are 
roughly at the same distance, the combination of shorter 
exposure time, higher Galactic absorbing column, and 
assumed lower gas temperature of the NGC 5170 halo, 
leads to less tight constraints on the hot halo properties 
for this galaxy. The lower temperature is particularly 
important, as Chandra is not very sensitive to gas at 
T < 0.2 keV. 

Although we have attempted to assume realistic values 
of T, Z, and rj for the halo gas, it is worth briefly illus- 
trating the impact of these assumptions on the derived 
quantities. In particular, tighter constraints would result 
for either galaxy if assuming a higher halo temperature or 
metallicity. For our best constrained case of NGC 5746, 
increasing both the assumed T and Z by 50% compared 
to the adopted values in Table [2] would reduce our upper 
limits on Lx and -Mhot by roughly 30%, while doubling 
the lower limit on (t coo \). For NGC 5170, the correspond- 
ing changes would be a 40% reduction of Lx and Mhot > 
with (tcooi) increasing by a factor 2.8. Also note that 
if employing uncertainties at the 90% confidence level 
rather than the more conservative 3tr limits in Table [H 
we would instead obtain Lx < 2.2 (4.1) x 10 39 erg s _1 and 
Mhot = 1.7 (2.5) xlO 9 M for NGC 5746 (NGC 5170). As 
will be apparent, none of these modifications would affect 
our overall conclusions. 
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4. COMPARISON TO SIMULATIONS OF GALAXY 
FORMATION 

The lack of any detectable extraplanar hot gas re- 
lated to stellar feedback is not surprising for either of 
our target galaxies, given the absence of any correspond- 
ing extraplanar Ha or radio emission. The derived con- 
straint on the hot halo luminosity of NGC 5746, however, 
is also an ord er of magn it ude b elow the naive expecta- 
tion from the iToft et all (|2002ft simulations for cosmo- 
logically accreted gas around a galaxy of this circular 
vel ocity. Care s hould be exercised in directly comparing 
the IToft et all (|2002D predictions to our results though, 
as the former also include the emission from within the 
optical disk. Moreover, as mentioned in § [H all the 
massiv e (v c > 250 km s _1 ) galaxies in the To ft et aTl 
(12002ft study were formed in the old "standard" CDM 
cosmology (il m = 1, Q\ = 0), with an assumed cos- 
mic baryon fraction /& = 0.05 or 0.10, well below cur- 
rent estimates (ft = 0.17; ISpergel et ai1l2007ft . and with 
no metal-dependence on the radiative cooling rate. The 
aim of this section is to explore the predictions of more 
sophisticated higher-resolution simulations assuming the 
currently prevailing ACDM cosmology, and test whether 
these can accommodate the observational constraints ob- 
tained for both our galaxies. 

4.1. Cosmological Simulations 

For this purpose, we have employed cosmological 
TreeSPH simulations of galaxy formation and evolution 
in a flat (fi m = 0.3, = 0.7) ACDM cosmology The 
cod e, described in detail in lSommer-Larsen et alj (2005) 
and I Romeo et al.l (|2006[ ), incorporates energetic stellar 
feedback in the form of starburst-driven gala ctic winds. 
Contr ary to our earlier simulations used by T oft et al.l 
(2002) to predict halo properties (as discussed in § [1] 
and § [2|) , the present simulations also incorporate self- 
consistent chemical evolution of the gas, including non- 
instantaneous recycling. They have also been evolved in 
time according to the entrop y equation solving scheme of 
Springel & Hernquist (2002), rather than being based on 
the thermal energy equation, a modification which pro- 
vides increased numerical accuracy in lower-resolution 
regions. In addition, the present simulations assume a 
cosmic baryon fraction of fb — 0.15 and have all been 
run from an initial redshift of z = 39 with a mass resolu- 
tion ei ght times higher than those described bv lToft et al.l 
(2002). At 0.5 Gyr before z = 0, each hot gas particle 
in the simulations was further split into 64 equal-mass 
particles to provide additional resolution of the hot gas 
phase, a procedu re very similar to the one described by 
iPeek et aTl (|2008h based on the same cosmological simu- 
lations. Th e final gas-phase resolution is thus 512 times 
that of the IToft et aTl (|2002ft simulations, with the best 
resolved individual galaxies at z = featuring a to- 
tal of ~ 1.3 x 10 5 hot gas particles within the central 
r = 100 kpc, a hot gas mass resolution of ~ 1.7 x 10 4 M , 
and an SP H smoothing lengt h of hspn « 1.5 kpc. As dis- 
cussed by IPeek et~a l. (2008), this resolution is sufficient 
to track thermal instabilities in the hot halo gas that may 
lead to the formation of warm (T ~ 10 4 K) clouds within 
the ambient halo. 

Although the simulation code includes stellar feedback, 
optionally enhanced to mimic AGN outflows, powerful 



starburst-driven winds only occur in the simulations at 
early times {z > 4-5), where they help to solve the angu- 
lar momentum problem, the missing satellites problem, 
and possibly other prob lems related to the cold dark mat- 
ter scenario (see, e.g., ISommer-Larsen et al.l l2003h . At 
z < 4, the SN II feedback is quite moderate in our sim- 
ulations; for a typical v c = 220 km s _1 disk galaxy, the 
star formation rate at z — is ^0.5 — 1 M Q yr _1 , com- 
parable to the modest star formation rate of the Milky 
Way disk. The feedback of the individual star particles in 
the simulations is furthermore spatially and temporally 
uncorrelated (this is particularly true for SN la, whose 
mechanical, but not chemical, feedback is consequently 
neglected), so gal actic winds do not e asily develop at low 
redshift (see also lDahlem et aTll2006l ). Hence, a compar- 
ison between the simulation output at z — and our two 
relatively quiescent galaxies seems justified. It is worth 
emphasizing that the feedback strength in the simula- 
tions is calibrated to reproduce the optical and morpho- 
logical properties of observed galaxies. The X-ray emis- 
sion is thus predicted without any additional adjustable 
parameters. 

For the calculation of X-ray properties of the simulated 
galaxies, we follo wed the general procedure outlined by 
IToft et aTl (|2002l ). For each of the simulated galaxies at 
z = 0, a catalog of SPH gas particle positions, densities, 
temperatures, and masses was created. Each SPH parti- 
cle was treated as an optically thin thermal plasma, and 
the associated X-ray luminosity calculated at the relevant 
position in a given phot on energy band us ing the meka 
plasma emissivity code (|Mewe et al.lll986ft . Total X-ray 
luminosities were then computed by summation over all 
particles in the volume of interest. The greatly increased 
numerical resolution of the present simulations at z = 
allowed us to bypass the additional s moothing over th e 
individual gas particles employed by IToft et aTl (|2002ft . 
However, as mentioned above, the increased resolution 
also facilitates the formation of small clouds within the 
hot halo via local thermal instabilities. Since the pres- 
ence of such relatively dense and cold clouds may un- 
physically boost the X-ray output of the hot gas in its 
immediate vicinity (due to the local smoothing of particle 
properties inherent in SPH techniques), hot gas particles 
within /isph of "cold" (T < 3 x 10 4 K) gas were discarded 
from the X-ray calculations. The removed gas generally 
constitutes an insignificant fraction (< 1%) of the total 
hot gas mass considered. 

A few comparisons between NGC 5746 and results from 
an earlier version of t hese simulations were performed in 
IPedersen et al.l (|2006l ). Here we present a more detailed 
comparison based on the updated Chandra analysis of 
NGC 5746, while also including the results for NGC 5170 
and for other, more actively star-forming, spirals as taken 
from the literature. Further comparison s are m ade to our 
previous results described in lToft et al.l |2002), with the 
X-ray properties of the latter here re-calculated for the 
adopted aperture and energy band. However, in order 
to provide as fair a comparison as possible between the 
two simulation samp les, we only consider the subset of 
the lToft et aTl (|2002ft galaxies that was formed formed in 
simulations assuming the same flat (f2 m = 0.3, J7a = 0.7) 
ACDM cosmology as adopted in our new simulations, 
and which assumed fi, = 0.10 (along with a fixed pri- 
mordial chemical gas composition). Note, as mentioned 
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Fig. 5. — 0.3—2 keV X-ray luminosities versus v c for our simu- 
lated galaxies (black squares), and for those of Toft et al. (2002; 
gray squares), all calculated assuming the NGC 5746 halo aper- 
ture. Also shown are observational res ults for extraplanar emission 
seen around other spira ls, taken fromlWa ng et all J2003I, "W03"), 
IStrickland et alj ll2004al "S04"1. ITullmann et al.l H2006al . "T06"), 
and lLi et aT . (2006, 2007, "L06+07"). 



above, that this excludes all the I Toft et alJ (|2002h galax- 
ies with v c > 250 km s , so a direct comparison of 
the Toft et al. results to our observed galaxies is not at- 
tempted here. 

4.2. Hot Halo X-ray Luminosities 

One of the hot halo properties that can be most 
straightforwardly compared to observations is the dif- 
fuse X-ray luminosity. For all simulated galaxies, we de- 
rived the halo Lx within the aperture also adopted for 
the NGC 5746 halo analysis (i.e. a spherical volume ex- 
tending to r = 40 kpc, with the L> 2 5 ellipse excluded as 
described in § 13. 3p . The results are presented in Fig- 
ure 03 which shows a clear trend in hot halo Lx with 
galaxy mass, measured by v c , for the simulated galax- 
ies. Also shown are the observational constraints on 
the halo luminosities of NGC 5170 and NGC 5746. The 
simulation results are generally consistent with the con- 
straints obtained for our observed galaxies, and we con- 
clude that, within the adopted apertures, the lack of de- 
tectable extraplanar X-ray emission around either target 
agrees with these updated predictions for emission from 
cosmologically accreted gas. 

There s eems to be reaso nable agreement between the 
results of iToft et all ((2002) and those of our newer sim- 
ulations, but with a tendency towards a somewhat lower 
Lx at a given v c for the latter. For example, for Milky 
Way-sized galaxies (v c in the range 200-250 km s _1 ), the 
mean halo luminosity for the IToft et ail (|2002f ) galaxies 
is 1.4 x 10 38 erg s _1 , 75% higher than the correspond- 
ing result for the more recent simulations. We attribute 
this in part to the increased numerical resolution of the 
latter, as discussed in more detail later. Nevertheless, 
the pseudo-bolometric (0.01-10 keV) halo luminosities 
Lx,boi in our new simulations remain consistent with 
the L x oc v%~ 7 relationship found by IToft et all (|2002D . 
with an orthogonal regression fit in logLx-logu c space 
yielding Lx,t>oi oc irj' 0±0 ' 7 for the new simulations. It 
is interesting to note that this dependence on the char- 



acteristic stellar velocity of the system mirrors that of 
the diffuse gas in early-type galaxies, whose X-ray lu- 
minosity roughly obeys Lx oc L? B cx <r 6 ~ 7 (based on 
the Faber- Jackson rel ation and the Lx~Lb relation of 
lO'Sullivan eTafll200lh . 

The amount of hot X-ray gas surrounding the galax- 
ies in our simulations depends not only on the mass of 
each galaxy's dark matter halo but also on the detailed 
accretion history of this halo. Numerical calculations in- 
dicate that accretion histories of galaxy-sized halos can 
be highly disp arate even for isolated galaxies of compara - 
ble mass (e.g. Ivan den Bosch 2002; D e Rossi et alT 2008). 
Hence, any dependence of hot halo properties on v c can 
a priori be expected to be subject to considerable scat- 
ter. Indeed, despite the clear trend in Lx with v c for 
the simulated galaxies, it is important also to empha- 
size the considerable variation in predicted halo Lx at a 
given circular velocity in Fig. 03 For example, the two 
most massive galaxies in our new simulation sample dis- 
play a factor of five variation in 0.3-2 keV Lx within 
the adopted aperture and yet are separated in "mass" 
by only Av c = 13 km s" 1 , comparable to typical obser- 
vational uncertainties on v c . This large variation at a 
given v c is clearly an issue to keep in mind when com- 
paring predicted and observed halo luminosities, or when 
attempting to predict hot halo fluxes for actual galaxies 
on the basis of simulation results. 

In order to also facilitate comparison of the simula- 
tion predictions to other observational results, we have 
compiled available constraints from the literature on ex- 
traplanar emission seen around other spirals, although 
we make no claims as to the completeness of this compi- 
la tion (note thoug h that we have not included the results 
of lSun et al.ll2007l because it is not entirely clear to what 
extent, if any, these spirals show extraplanar emission). 
These results are also included in Fig. 03 but we stress 
that the halo luminosities in these cases have not neces- 
sarily been derived within the same aperture as adopted 
for NGC 574 6 and our simula t ed gala xies. For example, 
the results oflStricklan d et ail (|2004al ) were all derived at 
a distance \z\ > 2 kpc from the disk midplane. If instead 
adopting this aperture, the increase in Lx predicted for 
our simulated galaxies is at the 10-40 per cent level, with 
no systematic dependence on v c . We expect similar cor- 
rections fo r the apertures used for the additional galax- 
ies in the ITullmann et al.l (|2006af) sample, where halo 
X-ray properties were constrained outside the Ha disk, 
again typically at \z\ > 2 kpc from the midplane. These 
changes are relatively small, well within the scatter seen 
for the simulated galaxies. 

The iLi et al.l (|2006l 120071 ) results in Fig. 02 include 
unresolved emission within regions of the optical disk 
(though with the prominent du st lane in the Sombrero 
galaxy excluded; ILi et al. 2007). Hence, the true halo 
luminosities of these galaxies should be somewhat lower 
than indicated in this figure, although the apertures em- 
ployed by Li et al. extend to smaller radii than adopted 
for NGC 5746 and hence our simulated galaxies. We 
include these two data points in Fig. 03 for complete- 
ness, in particular because the result for the massive 
(v c « 359 km s _1 ; HyperLeda database) Sombrero 
galaxy clearly poses a strong constraint on theoretical 
models, and since Fig. 03 suggests that the impact of its 
~ 10 9 Mq supermassive black hole on the diffuse Lx 
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must be rather modest. We also note that a very low 
halo X-ray luminosity of 3 x 10 38 erg s _1 has been re- 
ported forthe v c w 250 km s" 1 Sb spiral NGC4565 
(Wang 2005). This result would be broadly consistent 
with the new simulation predictions in Fig.[5]if assuming 
similar apertures. However, the result may be affected 
by emission fro m a nearby background cluster of galax- 
ies (|Wang 2005), and the details of the Chandra analysis 
have not been published, so we have chosen not to in- 
clude this data point in Fig. [5] 

Keeping the above caveats in mind, Fig. [5] indicates 
that the theoretical predictions can generally be accom- 
modated by the observational results for the (presum- 
ably mainly starburst-generated) hot halos seen in other 
studies. The latter show X-ray luminosities for low- and 
intermediate-mass spirals which are up to three orders 
of magnitude above the simulation results at a given v c 
but which are also roughly independent on v c and show 
considerably larger scatter at a given v c . It is, neverthe- 
less, i nteresting that X-ray observations of high-w c spirals 
(e.g., iLi et al.l l200l are now providing luminosity con- 
straints which are comparable to simulation predictions 
within the assumed aperture. 

4.3. Halo Gas Temperatures 

The inferred temperature of the hot halo gas in obser- 
vations and simulations is a useful diagnostic of the na- 
ture and origin of this gas. Figure [5] shows the emission- 
weighted temperature (T) of hot halo gas against v c , 
for all gas within the (100 kpc) 3 simulation volume sur- 
rounding each simulated galaxy. The predictions for the 
halo temperature are generally very similar to the ex- 
pected virial temperature T vlr cx v%, as anticipated for 
accreting, shock-heated material, but with a slight ten- 
dency for (T) to mildly exceed T V j r . A plausible expla- 
nation for this is that our estimate of (T) is dominated 
by emission close to the disk, where the gas temperature 
rises slightly above T vir , as demonstrated below. We note 
that a halo gas temperature slightly exceeding the virial 
temperature of the dark matter is also a generic feature 
of the phenomenolqgical h ot halo models considered by 
iFukugita fc Peebles! poM) . 

In contrast to the clear T-v c trend exhibited by 
the simulation s ample, published tem perat ures for the 
IStrickland etafl (|2004ah and lLi et all (|2006l . l2007t ) sam- 
ples show no obvious systematic dependence on v c 
(note that single emission-weighted halo tempe r atures 
ar e not available f o r the s amples of W ang et all (|2003l ) 
or lTullmann et all (|2006al ). so these galaxies are not in- 
cluded in this plot). Any trend with v c for the observed 
galaxies must thus be subject to considerable scatter, and 
is possibly modulated by a dependence on additional pa- 
rameters such as the star formation rate. Also note the 
absence of observed data points with T < 0.2 keV in the 
plot, which is at least in part related to the difficulty 
of detecting such low-temperature gas with Chandra or 
XMM. This also applies to the very low-w c M82 — ex- 
cluded he re for reasons of c l arity bu t included in Fig. [1] — 
for which [Strickland et al.l (|2004al ) list separate temper- 
atures of 0.27 and 0.37 keV for the halo regions covered 
by ACIS-S and ACIS-I observations, respectively. 

For the galaxies with detected extraplanar emission, 
the lack of a clear systematic trend in both Lx and T 
with v c , and hence galactic mass, strongly argues against 
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Fig. 6. — Global emission- weighted hot halo temperature versus 
v c for our simulated galaxies (squares), along with observational re- 
sults for extraplanar emission around other spirals, with symbols as 
in Fig. [5] Dashed line shows the virial estimate T v j r Ri (l/2)fim, p v^ 
of the hot halo temperature. 



a gravitational origin for this extraplanar emission. This 
interpretation gains support from the facts that the hot 
halo temperatures of these galaxies consistently exceed 
the virial temperatures expected from the circular ve- 
locity of the disk and are furthermore invariant with 
respect to Lx over four orders of magnitude in SFR 
(jGrimes et al.ll2005[ ). In general, the halo gas is almost 
certainly related to galactic wind or fountain activity in 
these cases, rather than representing externally accreted 
material. Indeed, unambiguous kinematic evidence for 
outflows perpendicular to t he disk has been uncov ered 
in many of these galaxies (|Strickland et al.l l2004bl and 
references therein). 

A single-temperature model was assumed when con- 
straining hot halo properties of our observed galaxies 
To test the validity of this assumption, we 



investigate spatial variations in the predicted hot halo 
temperature in Fig. [3 in the form of emission-weighted 
temperature profiles extracted perpendicular to the disk. 
The reasonably uniform temperature distribution pre- 
dicted outside the disk suggests that our assumption of 
a single-T model for the observed galaxies is well moti- 
vated. The one exception is the lowest-u c galaxy, which 
exhibits a very prominent warp in the cold gas distribu- 
tion that also affects the temperature distribution of its 
hot gas, giving rise to the abrupt change seen around 
\z\ « 25 kpc. We note that a qualitatively very similar 
picture to that in Fig. [7] applies to the metal abundance 
of hot halo gas, once again supporting the simplifying as- 
sumption of a single hot gas phase adopted in § I3.3I This 
is true irrespective of considering the emission- or mass- 
weighted abundances, although the latter is significantly 
lower within the adopted aperture, with a mean and stan- 
dard deviation of Z = 0.03 ± 0.01 Z Q compared to the 
corresponding emission-weighted value of 0.2 ± 0.1 Zq. 

Fig. [7] also shows that the predicted halo temperature 
is in good agreement with the virial temperature ex- 
pected for the simulated galaxies (i.e. T ss 0.3 keV for a 
galaxy with v c w 300 km s _1 ) at intermediate distances 
(\z\ ~ 20 kpc) from the disk, but drops slightly below 
this value further out. The overall halo X-ray emission 
is dominated by gas closer to the disk, where T gen- 
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Fig. 7. — Emission- weighted temperature profiles extracted per- 
pendicular to the disk for the simulated galaxies. The gray-scale 
reflects the expected virial temperature T v i r (darker colors indicate 
higher T v ; r ), with T v ; r itself marked for each profile by a filled cir- 
cle at an (arbitrary) large distance above the disk. Dashed vertical 
line marks the extent of the D25 minor axis of NGC 5746. 



erally increases gently to slightly above T vir , leading to 
the mildly super- virial global temperatures in Fig. [6] At 
even smaller distances from the disk midplane, the tem- 
perature shows evidence of a beginning decline as a sub- 
stantial fraction of halo gas eventually reaches average 
densities which allow for rapid cool-out. This is remi- 
niscent of a classical central cooling flow, but it is worth 
emphasizing that cool-out of halo gas in these simula- 
tions also proceeds via condensation of small pressure- 
supported clouds at various distances from the disk, with 
a monotonic increase in overall m ass cool-out rate with 
decreasing galact ocentric distance ([Sommer-Larsenl2006t 
iPeek et alll2008| ). Also note that although some gas is 
thus cooling rapidly in the central halo regions, Fig. [7] 
still suggests that the X-ray gas outside the disk could 
appear near-isothermal. This behavior could potentially 
compromise direct observational tests of whether halo 
gas is indeed cooling to below X-ray temperatures (e.g. 
via X-ray grating spectroscopy), as much of the cooling 
occurs very close to the optical disk where absorption 
by cold gas and contamination from other X-ray sources, 
diffuse and point-like, could be substantial. 

4.4. Hot Gas Mass and Cool-Out Rates 

In the context of understandin g the "missing baryon" 
problem in disk galaxies (e.g. [Mailer fc Bullock! 12004 
iFukugita fc Peebles! [20061 : ISommer-Larsenl[2006h . the to- 
tal amount of hot halo gas around such galaxies is a cru- 
cial quantity. For our simulations, the correlation be- 
tween hot gas mass Mhot within the NGC 5746 aperture 
and v c is displayed in Fig. [HJz. Given that the simulations 
also include cold, non-X-ray emitting gas around the 
galaxies and that Mhot is not an X-ray emission- weighted 
quantity, a low-temperature cut at T = 0.05 keV was 
imposed on the simulated galaxies in order to provide a 
sensible comparison to the Chandra constraints on X- 
ray emitting gas mass. This temperature cut is well 
below T V j r for all the simulated galaxies (cf. Fig. [6|), 
and is probably close to the lower limit at which gas 
can be reliably detected in emission in typical Chandra 
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Fig. 8.— (a) Mass of hot (T > 0.05 keV) func- 
tion of v c within the NGC 5746 halo aperture, with symbols as 
in Fig. \E\ (b) The corresponding mean hot gas cool-out rate, 

(M) SB Afgas/<*cool>. 



or XMM exposures (see , for example, the discussion in 
iRasmussen et al.ll2004bT) . 

As anticipated, a clear trend with v c is present for 
the simulation results. These are also seen to be con- 
sistent with our observational constraints, but we note 
that the 3<r upper limit derived for NGC 5746 is here 
only within a factor of 2-3 of the simulation predictions 
for a galaxy of this mass. Clearly, any changes in the 
assumed input physics in the simulations that would re- 
sult in significantly higher hot halo masses at z = 
are disfavored by our observations. Also note that the 
overall hot gas masses within this aperture at a given 
v c > 200 km s" 1 are similar for the two simulation sam- 
ples, despite the clear differences in corresponding Lx 
visible in Fig. [5] The main explanation is that the in- 
creased numerical resolution employed here, and the re- 
moval of a small amount of hot gas particles within an 
SPH smoothing length of cold gas in the new simula- 
tions, reduces the amount of (artificially) strongly X-ray 
emitting high-density halo gas. 

The predicted amount of hot halo gas within the 
adopted aperture is Af ho t < 10 9 M Q for Milky Way- 
sized galaxies, corresponding to only < 2% of the bary- 
onic disk mass of such galaxies. The situation for our 
observe d galaxies is unlikely to be significantly different. 
In fact, iRand fc BeniamirJ (|2008f ) used the inferred Hi 
rotation curve of NGC 5746 to derive a dynamical mass 
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inside r « 40 kpc of 7.9 x 10 11 M Q , which, combined with 
Mhot from Table [21 implies that any hot halo gas outside 
the optical disk contributes less than 0.3% of the total 
mass budget within this radius. In terms of baryonic 
mass, the inferred Hi mass in the disk of 1 x 10 10 M , 
added to a stellar mass of 1.6 x 10 11 M Q as estimated 
from the o ptical parameters in Ta ble Q] using the pre- 
scription of iMannucci et all (|2005l ). implies that < 1% 
of the baryons in NGC 5746 are in a hot phase within 
r = 40 kpc. For NGC 5170, the corresponding constraint 
is < 4%, based on a stellar mass of 6 x 10 10 Mp and an 
Hi mass of 1.4 x 10 10 M Q (jBottema et al.lll987l ). 

These fractions are expected to grow drastically with 
galactocentric distance, however, as most of the diffuse 
hot baryons associated with our simulated galaxies re- 
side beyond the r = 40 kpc aperture probe d here. As 
discussed in detail in ISommer-Larsenl 12006;. the total 
mass of hot halo gas within the virial radius may equal 
or exceed the mass of cold baryons in the disk of a simu- 
lated v c w 300 km s _1 galaxy. For our observed galaxies, 
this possibility is still allowed by our constraints in Fig- 
ure [H provided that the radial distribution of any hot 
halo gas around our galaxies is not significantly steeper 
than in these simulations. Deeper X-ray observations 
would be required to test this scenario in detail, but ob- 
servational verification of this prediction would clearly 
help in resolving the well-known "missing baryon" prob- 
lem in disk galaxies. 

As a final comparison between the simulations and our 
observational constraints, Fig. [8j> shows the results for 
the mean cool-out rate (M) m Mhot/(*cool) of hot gas 
within the NGC 5746 aperture. Again, a temperature 
cut at T — 0.05 keV was imposed on the simulations, 
but M was otherwise estimated in the same manner as 
for the observed galaxies (§ 13 .3(1 . i.e. with hot gas cooling 
times (tcooi) derived a ssuming the metal-dep endent cool- 
ing curves A(T, Z) of ISutherland fc Dopital (|1993t ). The 
figure reveals substantial scatter at any given v c and no 
clear systematic variation with this quantity. The sim- 
ulation predictions again remain consistent with our ob- 
servational r esults as li s ted in Table [2] The difference 
between the iToft et all (|2002j ) and newer simulations is 
largely due to lower mean cooling times for the hot gas 
in the former, for the reason discussed above. 

It is important to note, however, that the estimated 
values of M should not be identified with the actual ac- 
cretion rate of hot halo gas onto the disk, in part be- 
cause these estimates do not include the gas condens- 
ing out of the hot halo p hase due to small -scale ther- 
mal instabilities (see, e.g. iPeek et al"1l2008l) . but more 
importantly because our M does not take into account 
the gas cooling inside the inner boundary of the adopted 
aperture. In practice, the derived values of M should 
be considered lower limits to the true disk accretion 
rates in the simulations. For a detailed discussion of ac- 
tu al accretion rates for the sim ulate d galaxies, we refe r 
to ISommer-Larsen et al.l (|2003t ) and IPeek et al.l (|2008h . 
Such a discussion is beyond the scope of this work, also 
because it is not straightforward to compare these rates 
to our observational constraints in Table [2j 



5. DISCUSSION 



We emphasize that the observational part o f this study 
was in itiated on the basis of the predictions of lToft et al.l 
( 2002) , which suggested that the performed observations 
of NGC 5746 would be sufficiently sensitive to reveal the 
presence of an accreted hot halo. Our newer simula- 
tions predict somewhat lower hot halo luminosities out- 
side the disk, in consistency with the inferred absence 
of detectable hot halo gas around both NGC 5746 and 
NGC 5170. It is nevertheless clear that observational re- 
sults for massive galaxies are within factors of a few of 
the theoretical predictions for halo luminosity and gas 
mass, and can potentially begin to test the input physics 
in cosmological simulations. For example, the predicted 
hot gas mass for v c rj 300 km s _1 galaxies is within a fac- 
tor of 2-3 of the upper limit derived for NGC 5746. Nu- 
merical effects aside, inferred halo X-ray luminosities of 
even more massive galaxies such as NGC 4594 (cf. Fig. [3]) 
clearly leave limited scope for changes in the simulation 
input physics that would raise the predicted present-day 
halo Lx, though with the caveat that the simulations 
show some "cosmic variance" at a given v c . Further high- 
sensitivity observations may begin to seriously challenge 
the theoretical models and provide useful constraints on, 
for example, the galactic feedback prescriptions imple- 
mented in state-of-the-art numerical work. 

Com pared to the earlier results presented bv lToft et al.l 
(|2002j ). it is also clear that the X -ray luminosity and cool- 
out rate of hot halo gas are somewhat lower in the more 
recent simulations. One possible reason is the higher 
mass resolution of the latter, which helps to suppress 
the importance of a problem associated with the spa- 
tial smoothing inherent in SPH techniques, namely the 
possibility that the properties of hot gas particles can 
become strongly affected by those of much colder and 
denser particles in their immediate vicinity. As a rough 
test of the impact of improving the hot gas mass resolu- 
tion, we re-simulated our lowest-mass (v r . ~ 180 km s -1 ) 
galaxy from z = 39 at 64 times the IToft et al.l (|2002l ) 
resolution, with each hot gas particle then subsequently 
split into eight equal-mass ones at a later stage. At high 
z, this version is thus eight times better resolved than the 
one shown in Fig. [5j and the halo Lx is here found to be 
50% lower. An additional test for a more massive galaxy, 
the z w 1 version of the v c « 245 km s _1 disk in Fig. [5l 
showed little change in halo Lx at higher resolution, how- 
ever. This suggests that higher resolution may decrease 
the predicted Lx somewhat, but the reduction is moder- 
ate (well within the factor 4-5 in "cosmic variance" that 
is evident in Fig. [5]), and Lx may be close to conver- 
gence at the high-mass end. We do not currently have 
higher-resolution versions of our most massive galaxies 
at z = due to the prohibitive CPU times required, so 
more detailed tests of the impact of numerical resolution 
will remain the subject of future work. 

Two other effects also act to reduce the present-day 
hot halo luminosity in the updated simulations, namely 
the higher assumed cosmic baryon fracti on fb = 0.15 (vs. 
fb = 0.10 for the comparison galaxies of lToft et al.ll2002l 
considered here) and the higher metal abundance result- 
ing from implementing chemical evolution in the simu- 
lation code. This may at first seem counter-intuitive, 
but the increased radiative cooling efficiency implied by 
both these modifications actually results in less hot and 
dense gas remaining in the halo at z — for a given v c 
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(see also iRasmussen et al.ll2004atlSommer-Larserill2006l) . 
Note that both the hot gas fraction and the total "exter- 
nal" (non-disk) baryon fraction inside the virial radius 
are reduced at fixed v c when assuming higher /;,, not 
only because the accreted shock-heated gas can leave the 
X-ray phase more rapidly, but also because the virial 
radius itself decreases at a given v c as a consequence 
(|Sommer-Larsenll2006D . 

Despite the reduced X-ray output in the current sim- 
ulations, we readily ack nowledge the points raised by 
iLi et all (|2006l I2007T ) and I Wand (|2007h . Their compari- 
son of observational results obtained for NG C 2613 and 
NGC 4 594 to the 0.2-2 keV luminosities of iToft et all 
(2002) indicate that the latter overpredict total (disk 
plus halo) intrinsic X-ray luminosities by at least a 
factor of 10 for these massive galaxies. This discrep- 
ancy is at least in part related to the common over- 
cooling problem in hydrodynamical simulations, a prob- 
lem w hich also plagues t he current simulations at some 
level (|Romeo et al.l l2006l ). However, driven by the re- 
quirement to avoid confusing any halo emission with 
emission produced by unrelated processes, we have to 
some extent circumvented that issue here by only consid- 
ering the X-ray propertie s of gas out si de th e optical disk. 
Note also that both the IToft et al] (|2002ft and present 
simulations do not incorporate the effects of intrinsic X- 
ray absorption by cold gas and dust in the disk, which can 
naturally be substantial for edge-on galaxies and may be 
non-trivial to properly correct for across the full optical 
extent of the galaxy in observational studies. 

While our results suggest that increased numerical res- 
olution may help to alleviate the issue of ove r-cooling 
somewhat, i t is also possible, as hypothesized bv lLi et all 
( 20061 l2"007h . that some over-cooling in the simulations is 
related to incomplete treatment of feedback from SN la, 
if these can provide large-scale distributed heating. As 
mentioned, mechanical feedback fro m SN l a is no t in- 
cluded in the present simulations, but IWand (|2007l ) sug- 
gest that numerous SN la in early-type spirals (such 
as NGC 2613, 4594, and 5746) could generate a bulge 
wind which would be sufficiently powerful to suppress 
the cooling of any infalling halo material and lower its 
X-ray luminosity. Observational verification of the pres- 
ence of SN la-driven winds in spirals would certainly 
provide an important new ingredient to be incorporated 
in numerical simulations. Our X-ray data of NGC 5746 
and NGC 5170 do not allow a direct observational test of 
this scenario, apart from requiring such outflows to have 
X-ray luminosities below the inferred halo upper limits. 
The discussion pertaining to Fig. [IJ although based on 
somewhat simplistic assumptions, suggests that direct 
supernova blow-out of disk material is unlikely to be im- 
portant in our two galaxies even if including SN la. Our 
galaxies may therefore not represent optimal targets for 
a n observation al test of this scenario. 

IWand (|2007f ) also proposes another possible solution to 
the over-cooling problem, noting that the halo gas may 
be chemically inhomogeneous but generally very low in 
metallicity. While metal-rich clouds can cool out rapidly 
and fall towards the disk, most of the halo gas may thus 
radiate inefficiently and would be difficult to detect in 
either emission or absorption. Although an interesting 
possibility, it is not obvious that this explanation may 
necessarily help to solve the over-cooling problem. As 



discus sed above, the metal-free halos in the IToft et ahl 
(2002) simulations are actually more X-ray luminous at 
a given v c at zero redshift than those in our more recent 
simulations. Other factors such as the assumed cosmic 
fb in the simulations may play a role for this result, as 
already mentioned, but this difference between metal- 
enriched and metal-free halos persists even when consid- 
ering halos in simulations run with iden tical cosmological 
param eters and numerical resolution ([Rasmussen et aTl 
l2004al) . Allowing for chemical enrichment in our simu- 
lations systematically reduces the predicted present-day 
halo luminosity for a given v c . 

On the observational front, it is generally well es- 
tablished that disk galaxies have experie nced consider- 
able stellar mass growth since z ~ 1 (e.g. iBarden et al.l 
120051 and references therein), and both observational 
and theoretical evidence suggests that gas consumption, 
rather than major merger activity, has played the dom- 
inant role in d riving this growth (Ide Ravel et all 120081 : 
iDe Rossi et~al1 120081 : iDaddi et al.ll2008h . The hot halos 
discussed in this paper could potentially have supplied 
most of the fuel required for the continued growth of 
disks, provided that cool-out of such gas is not strongly 
inhibited by, for example, disk feedback processes. De- 
spite the absence of a detectable large-scale X-ray corona 
around either of our target galaxies, hot coronal gas sur- 
rounding both targets could still be present below our 
detection limits. It remains unclear, however, to what 
extent any such gas is currently cooling out to replenish 
the gas consumed by star formation in their disks, and 
if so, whether this proceeds mainly as a central cooling 
flow or through condensation and infall of small clouds 
across a wide range of galactocentric distances within a 
hot halo. For NGC 5746, our naive constraint on the 
cool-out rate of hot halo gas outside the optical disk, 
M < 0.3 M Q yr -1 , suggests that large-scale cooling 
from a hot halo cannot balance the SFR in the disk of 
0.9 ±0.2 M© yr -1 . It is important to reiterate, however, 
that the estimated cooling rates should be considered 
lower limits to the total rate at which halo gas may leave 
the hot phase and accrete onto the disk (§ 14. 4p . so our 
results are not necessarily inconsistent with the possi- 
bility that a hot halo around NGC 5746 can supply the 
disk with material at rates comparable to the SFR. For 
NGC 5170, the weaker upper limit on M is a factor of ~ 3 
above the inferred SFR, again precluding any definitive 
conclusions in this regard. 

As suggested by a number of theoretical models, and, 
indee d, by the simulations presented here (see lPeek et alJ 
2008i ). local thermal instabilities within a hot halo can 
also produce small pressure-supported clouds that would 
fall towards the disk on time-scales well below the cool- 
ing time of the ambient hot gas. This would contribute 
to feeding ongoing star formation and could potentially 
represent the dominant accretion mo de for Milky Way- 
sized spirals ([Mailer fc Bullo ck 2004). However, Hi ob- 
servations of nearby spirals generally suggest that infall 
of cold gas does not occur at rates sufficient t o sustain 
star formation in the disk ([Sancisi et all 120081 ). Hence, 
infall of additional, possibly hot, material (at an aver- 
age rate of « 1 M yr -1 ) seems required to balance gas 
co nsumption by sta r formation, in line with the findings 
of iPeek et al.l (|2008h . Of particular relevance here is the 
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study of lRand fc Beniaminl (|2008l ). who searched for ev- 
idence of an extended Hi component around NGC5746 
associated with infalling cold clouds. While the presence 
of vertically extended Hi was detected, the majority of 
this gas is best interpreted as a warp rather than as an ex- 
tended Hi halo such as has b een reported around a num - 
ber of other nearby spirals ([Fraternali fc Binnevl [2006) . 
A number of high-|z| neutral clouds were also seen, pos- 
sibly representing infalling material, though association 
with a disk-halo flow cannot be excluded for most of 
these. No corresponding study has, to our knowledge, 
been conducted for NGC5170. 

Thus, given the present data, it remains debatable 
whether the current SFR in galaxies such as NGC 5746 
can be balanced by infall of halo material. Cold gas ac- 
cretion alone seems insufficient both in NGC 5746 and in 
other nearby spirals, and our upper limit on the hot gas 
infall rate does not include gas cooling out within the 
disk region itself, so the interpretation of the observa- 
tional X-ray results is not straightforward. It is apparent, 
however, that cooling of hot halo gas outside the optical 
D25 ellipse may not be abl e to balance star formation 
in NGC 5746. Interestingly, iFraternah fc Binnevl (|2008l ) 
offer a different accretion scenario that may reconcile the 
discrepant accretion and SFR estimates in nearby spirals 
without invoking substantial direct cool-out of hot halo 
material. In their model, accretion has an important con- 
tribution from supernova-ejected fountain clouds which 
sweep up (hot or cold) ambient halo gas before return- 
ing to the disk, losing angular momentum in the process. 
Their combined fountain plus accretion models provide 
a better match to the extraplanar Hi velocity distribu- 
tion in a few well-studied edge-on spirals than a fountain 
model alone, and yield accretion rates comparable to the 
SFR in these galaxies. Direct observational verification 
of this scenario may prove difficult, however, and cosmo- 
logical simulations still lack the resolution to adequately 
capture such processes. 

6. SUMMARY AND CONCLUSIONS 

We have used X-ray data to re- visit our earlier claimed 
detection (|Pedersen et al.l l2006f) of a large-scale dif- 
fuse X-ray halo surrounding the massive edge-on spiral 
NGC 5746. An analogous analysis of the less massive 
NGC 5170 has been completed for comparison. Our up- 
dated Chandra analysis employs calibration data which 
accounts for the positional dependence of the contami- 
nant on the Chandra optical blocking filters (unlike the 
iPedersen et al.l [20061 analysis, which could not take this 
into account), and has failed to reveal a statistically sig- 
nificant detection of extraplanar X-ray emission around 
either galaxy. We find that both the spatial and spec- 
tral hot halo p roper ties of NGC 5746 as reported by 
IPedersen et al.l (|2006l ) are consistent with being due to 
the spatial variation of the degradation in the ACIS CCD 
quantum efficiencies. 

Both our target galaxies show no signs of significant 
nuclear or star formation activity. This is supported 
by our Ha imaging of both galaxies and by radio data 
taken from the literature, suggesting that the constraints 
on extraplanar X-ray emission obtained for both galax- 
ies can be used to also place constraints on the prop- 
erties of intergalactic, shock-heated material accreting 
onto their disks as predicted by disk galaxy formation 



models and cosmological hydro-simulations. If assum- 
ing halo gas close to the expected virial temperature as 
suggested by such theoretical work, we constrain the to- 
tal 0.3-2 keV luminosity of the halos within r — 40 kpc 
from the optical disk to L x < 4.0 x 10 39 (NGC 5746) and 

< 7.5 x 10 39 erg s" 1 (NGC 5170) at 3a significance. The 
corresponding hot gas masses are Mhot < 2.3 x 10 9 and 

< 3.4 x 10 9 M Q , respectively. 

We have performed a detailed comparison of these 
constraints to results of recent cosmological simulations 
of galaxy formation and evolution. These s imulations 
repres ent an update of our earlier work in iToft et all 
(2002), and predict hot halo luminosities and gas masses 
which are consistent with the limits obtained for both 
NGC 5746 and NGC 5170. Outside the optical disk, the 
simulations show h ot halo luminosit ies which are some- 
what lower than the lToft et a l. (2002) results, specifically 
by a factor of two for Milky Way-sized galaxies, mainly 
due to increased numerical resolution, the assumption of 
a higher baryon fraction, and the inclusion of chemical 
feedback. We find that "cosmic variance" in simulated 
halo properties can be quite substantial, with halo Lx 
potentially varying by a factor 4-5 for galaxies of com- 
parable mass as measured by their disk circular velocity 

Re- 
published results for all other galaxies with detected 
extraplanar emission can generally accommodate our 
simulation predictions, but observational results for high- 
mass spirals are clearly providing strong constraints on 
theoretical models for disk galaxy evolution. Based on 
the inferred halo temperatures and luminosities, we con- 
firm that the hot halo emission so far observed around 
other spirals is in general predominantly related to ac- 
tivity in the disk rather than being caused by infalling, 
externally accreted material. 

An unambiguous detection of externally accreted X- 
ray emitting gas around spirals is still pending. Deeper 
X-ray observations of carefully selected targets are re- 
quired to detect such halos, but, encouragingly, our sim- 
ulations suggest that this remains within reach of current 
X-ray instrumentation. Based on a rather modest Chan- 
dra exposure of less than 40 ks, we have constrained an 
X-ray luminosity and hot gas mass around NGC 5746 
which arc both within a factor of a few of the simulation 
predictions for a galaxy of t his mass. Oth er existing ob- 
servational constraints (e.g-. lLi et al.ll2007l ) are also close 
to the theoretical results and can potentially begin to 
test the assumed input physics in numerical simulations 
of galaxy formation and evolution. 
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